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ABSTRACT 

We use spatially resolved long-slit spectroscopy from Magellan to investigate the extent, kinematics, and ion- 
ization structure in the narrow-line regions of 15 luminous, obscured quasars with z < 0.5. Increasing the dynamic 
range in luminosity by an order of magnitude, as well as improving the depth of existing observations by a sim- 
ilar factor, we revisit relations between narrow-line region size and the luminosity and linewidth of the narrow 
emission lines. We find a slope of 0.22 ± 0.04 for the power-law relationship between size and luminosity, sug- 
gesting that the nebulae are limited by availability of gas to ionize at these luminosities. In fact, we find that the 
active galactic nucleus is effectively ionizing the interstellar medium over the full extent of the host galaxy. Broad 
300- 1000 km s"') linewidths across the galaxies reveal that the gas is kinematically disturbed. Furthermore, 
the rotation curves and velocity dispersions of the ionized gas remain constant out to large distances, in striking 
contrast to normal and starburst galaxies. We argue that the gas in the entire host galaxy is significantly disturbed 
by the central active galactic nucleus. While only ~ 10^-10^ Mq worth of gas are directly observed to be leaving 
the host galaxies at or above their escape velocities, these estimates are likely lower limits because of the biases in 
both mass and outflow velocity measurements and may in fact be in accord with expectations of recent feedback 
models. Additionally, we report the discovery of two dual obscured quasars, one of which is blowing a large-scale 
10 kpc) bubble of ionized gas into the intergalactic medium. 



1. INTRODUCTION 

It has become fashionable to invoke feedback from accreting 
black holes (BHs) as an influential element of galaxy evolution 
(e.g iTabor & Binnevlll993l; ISilk & Reeslll998l: iHopkins et al.1 
l2ooa ls ironi & SocratesI |20I0|) . Regulatory mechanisms are 
sorely needed to keep massive galaxies from forming too many 
stars and becoming over l y massive or blue at late times (e.g., 
iThoul & Weinber j1l995UCroton et al.ll2()06i) . Feedback from 
an accreting BH provides a tidy solution. For one thing, the 
gravitational binding energy of a supermassive BH is com- 
pletely adequate to unbind leftover gas in the surrounding 
galaxy. Furthermore, using simple prescriptions for black hole 
feedback leads to a natural explanation for the observed scal- 
ing relations between the BH mass and properties of the sur- 
rounding galaxy, including stellar velocity dispersion and bulge 
luminosity and mass (e.g., iHaring & Rixll2004 '; Giiltek in et al.l 
12009) . The problem remains to find concrete evidence of BH 
self-regulation, and to determine whether or not accretion en- 
ergy has a direct impact on the surrounding galaxy. 

There are some special circumstances in which accretion en- 
ergy clearly has had an impact on its environment. For in- 
stance, jet activity in massive elliptical galaxies and bright- 
est cluster galaxies d eposits energy into the hot gas enve- 
lope (see review in iMcNamara & NulsenI |2007|) . although 
the efficiency of coupli ng the accretion energy to the gas 
remains uncertain (e.g., iBirzan et all l2004i iBest et all 120051; 
iHeinz et al] l2006l) . as does the relative importance of heat- 
ing by the active nucleus as opposed to other possible 
sources (e.g.. IZakamska & NaravanI [20031: ISiiacki et al.l 120081: 



IConrov & OstrikeH 120081 : iParrish et al.]|2009l ). Likewise, there 
is clear evidence that powerful radio jets entrain warm gas and 
carry significant amounts of material ou t of th eir host galaxies 
(e.g., [v an Breugel et al. 1986; Tadhunte ^[T99Tt [W hittle 19923 
IVillar-Martm et al.l 1 19991: jNesvadba et al. 2006, 2008). How- 
ever, as only a minority 10%) of active BHs are radio-loud, 
invoking this mechanism as the primary mode of BH feedback 
would require all galaxies to have undergone a radio-loud phase 
- a conjecture which lacks direct evidence and contradicts a the- 
oretical paradigm in which r adio-loudness is determine d by the 
spin of the black hole (e.g.. lTchekhovskov et al.ll20I0l) . Thus, 
it is not clear whether BH activity in radio galaxies a ccounts 
for more than a small fraction of the BH growth (e.g., ISoltanI 
1982; Yu & Tremaine 2002) and therefore whether this mode 
of feedback is in fact the dominant one. 

Nuclear activity is known to drive outflows on small scales. 
Broad absorption-lin e troughs are seen in ~ 10-20% of lu- 
minous quasars (e.g. iReichard et al]|2003l) . and there is good 
reason to believe that the outflows are ubiquitous but have a 
cover ing fraction of ^ 20%, at least for high-Lboi/LEdd systems 
(e.g.. I Wevmann et alj [19811 [Gallagher et al.1 120071; [Shen et al.l 
2008). The velocities in broad absorption lines are high (~ 
10, 000 km s"' ), and they most likely emerge from a wind blown 
off of the accretion disk (e.g., Proga & Kallman 2004) . In a few 
rare objects the out flow appears to e xtend out to large distances 
from the nucleus (Moe et ani2009f) . but it is unclear whether 
most of these outflows have any impact beyond hundreds of 
Schwarzschild radii. Narrow associated absorption-line sys- 
tems are signposts of outflows extending to larger distances, 
but determining 
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Table 1. Sample and Observations 



SDSS Name 




log i|oiii| 


Obs. Date 




N 


1, 


PA 




F] 4 GH; 


(1) 


(2) 


(3) 


(4) 




|5| 


(6| 


(71 


(8) 


{9) 


(!0| 


SDSS JO 157 16.92-005304.6 


0.422 


42.77 


2007-08- 


15 


4500 


1 


180 


88 


75 


<0.92 


SDSS J084 135.09+0 10 156.3 


0.111 


42.44 


2007-03- 


IS 


7200 


2 


60 


44,74 


1 


3.99 


SDSS Jl 10621. 96-M335747.1 


0.242 


42.59 


2007-03- 


18 


5400 


2 


60 


60,149 


58 


<0.97 


SDSS Jl 12437.40+045618.8 


0.283 


42.66 


2007-03- 


19 


3600 


1 




37 


--90 


918.4 


SDSS Jl 14258.39+102747.8 


0.223 


42.81 


2007-03- 


19 


3600 


1 




119 




1.50 


SDSS J 1222 17.85-000743 .7 


0.173 


42.88 


2007-03- 


18 


7200 


2 


60 


84.174 


103 


4.24 


SDSS J 125334.49-034 158.1 


0.239 


42.42 


2007-08- 


14 


1800 


1 


120 





-37 


2.68 


SDSS J135646.10+102609.0 


0.123 


42.77 


2007-03- 


18 


3600 


2 


60 


179,134 




59.6 


SDSS J141315.30-014220.9 


0.380 


42.74 


2007-08- 


14 


2700 


1 


600 





46 


3.07 


SDSS J150719.94+O029O5.O 


0.182 


42.56 


2007-03- 


18 


4500 


1 


60 


29 


71 


4.55 


SDSS J205629.85+O05752.6 


0.435 


42.71 


2007-08- 


15 


5400 


2 


600 


109, 71 


182 


<2.4 


SDSS J2 12653.82+003553.0 


0.235 


42.27 


2007-08- 


14 


4788 


2 


120 


0,90 


-78 


9.3 


SDSS J22 12 17.92-094407. 8 


0.181 


42.26 


2007-08- 


14 


1800 


1 


120 


159 




<2.4 


SDSS J224027.05+004347.4 


0.438 


43.02 


2007-08- 


15 


3600 


1 


600 


129 


166 


<0.75 


SDSS J233558.39-OO5O0O.1 


0.439 


42.78 


2007-08- 


15 


3600 


1 


600 


134 


179 


3.99 



Note. — Col. ( 1): Official SDSS name. Col. (2): Redshift. Col, (3): L[oiii| (erg s^') as measiired from the SDSS spectrum of 
Reyes et al. 2008; used in the object selection. Col. (4): Date of LDSS3 observation. Col. (5): Total duration of spectroscopic 
exposure (sec). Col. (6): Number of slit positions. Note there are a couple of spectra from Paper I that were not useful for our 
purposes here, and thus are excluded. Col. (7): Exposure time (sec) for acquisition images. Col. (8): Position angle (deg) of 
each slit position, measured E from N. Col. (9): Posidon angle (deg) of galaxy, based on our fits in Paper I. PA is listed for the 
most extended galaxy component m all eases (e.g., the disk, if there is one). In the most disturbed cases, no PA is listed. Col. 
(10): Flux density (mJy) at 1.4 GHz measured either fn>m the FIRST survey or from the NVSS survey in cases when FIRST 
coverage is not available (SDSS J2056+0057 and SDSS J2212-0944) or the source is faint and extended (SDSS J2126+O035). 
SDSS Jl 124+0456 (4C +05.50) is a double-lobed radio galaxy (position angle 141 degEof N). 
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Table 2. Nebular Properties 



SDSS Name 


^'-'6 ^[OIII] 


lug i\[OIII] 


lug iVcont 


lug 1 v¥ mvifoill] 


R 










f5^ 




JOl 57-0053 


42.91 


<4.3 


4.0 


3.0 


2.9-3.2 


J0841+0101 


41.88 


4.0* 


4.0 


2.7 




Jl 106+0357 


42.52 


4.0 


4.1 


2.8 


2.9-3.2 


Jl 124+0456 


42.03 


4.3 


4.5 


2.6 


1.6-3.4 


Jl 142+1027 


42.49 


3.6* 


4.2 


2.6 




J 1222-0007 


42.63 


4.1* 


3.8 


2.7 


2.0-3.3 


J 1253-0341 


42.31 


4.2 


4.0 


2.6 


2.2-2.6 


J1356+1026 


42.23 


4.0* 


4.0 


2.6 




J1413-0142 


42.25 


3.8 


3.8 


2.7 


3.7-3.8 


J 1507+0029 


42.11 


3.8 


4.2 


2.8 


2.3-2.8 


J2056+0057 


42.60 


4.3 


4.0 


2.9 


2.7-3.2 


J2 126+0035 


42.06 


3.9 


4.1 


2.7 


2.3-3.2 


J22 12-0944 


41.39 


4.2 


4.0 


2.7 


2.5-2.6 


J2240+0043 


43.22 


<4.1 


3.9 


2.9 


2.9-3.3 


J2335-0050 


42.85 


<4.3 


3.8 


2.8 


3.2-3.3 



their physical radii (and thus the mass outflow rate) is notori- 
ously challenging. In cases where it is possible, the estimated 
outflow rates are thought to be sig nificant fractions of the ac- 
cretion onto the BH (see review in Crenshaw et al. 

It is clear that some quasars affect their environment some 
of the time. The extent and the dominant mode of these in- 
teractions remain open to interpretation. In particular, it is not 
clear whether quasars are effectively removing the interstellar 
medium (ISM) of their host galaxies during the high accretion 
rate episodes - those that account for the majority of the BH 
growth. Such feedback has been p ostulated by numerical sim- 
ulations (e.g., ISpringel et an2005h , but direct observational ev- 
idence for this process is lacking. 

In this work, we look for direct evidence of extended warm 
gas in emission, using the narrow-line region (NLR) and specif- 
ically the strong and ubiquitous [O III] A5007 line. The 
NLR is in some respects the ideal tracer of the interface 
between the galaxy and the active galactic nucleus (AGN), 
as the gas is excited by the AGN but extended on galaxy- 
wi de scales. For a long time, following the seminal work 
of IStockton & MacKermT ( Il987h , it was thought that truly ex- 
tended emission-line regions (so-called EELRs with radii of 10- 
50 kpc) were only found in radio-loud objects. Using narrow- 
band imaging, these authors examined known luminous, z ~ 0.5 
AGNs and found that ^ 1/3 of the radio-loud objects had lumi- 
nous extended [O III] nebulosities, while none of the radio-quiet 
objects did. It is not clear if the extended gas has an internal or 
external origin nor whether it is only present in radio-loud sys- 
tems or is only well-illuminated in the presenc e of radio jets 
fe.g., IStockton et al.ll2006l:lFu & Stocktonll2009h . 

Emission-line regi ons around radio-quiet systems 
dHusemann et alj|2008 f) are not usually as extended nor as lu- 
minous as those seen in the presence of powerful radio jets. 
This statement depends on the flux limit. At very low surface- 
brightness levels (^ 10~"*erg s"' cm"' arcsec"^), diverse mor- 
phologies are observed in emission line gas (e.g., Colbert et al. 
Il996t IVeilleux et al.l l2003). An interesting exception may be 
the broad-line active galaxy Mrk 23 1 . This galaxy shows out- 
flowing neutral and ionized gas that is exte nded on ~ 10 kpc 
scales and moving at thousands of km s"' (iHamilton & Keell 
[19871) . There is a jet in this galaxy (as weU as a starburst) but 
the jet i s not likely the so urce of acceleration of the neutral 
outflow dRupke et al.ll2005h . 



Note. — Col. (1): SDSS name. Col. (2): Nuclear L[oiii] (erg s"'; within 2'.'25) . 
Col. (3): Radius of [O ill] nebulosity (pc), as measured from the long-slit spectrum (§3). 
Upper limits are labeled, and * indicates double continuum sources, for which the size 
measurement may not be completely reliable. Col. (4): Radius of galaxy continuum 
(pc), measured in the same way as the [O in] sizes (§3). Col. (5): Full-width at half- 
maximum (FWHM) of the nuclear [O ill] spectrum (km s~' ). Col. (6): Power-law drop 
of luminosity in the outer parts of the nebula; L[oni] . Range includes both sides 

of all spectra for a given object. 

Rather than focus on unobscured (broad-line) quasars, where 
detailed study of the NLR extent and kinematics is hampered 
by the presence of a luminous nucleus, we look instead at ob- 
scured quasars. The experiment is worth revisiting in light 
of the discovery of a large sample of obscured quasar s with 
the Sloan Digital Sky Survey (SDSS; lYork et aLlBOOOl) . The 
sample, wit h z < 0.8, was select ed based on the [O III] line 
luminosity (Zakamska et al."2003) and now comprises nearly 
1000 objects (Reyes et al . 2008). E xtensive foll ow-up with the 
Hubble Space Telescope (Zakams ka etal.l 20061), Cha ndra and 
XMM-Newton ('Ptak et al. 200^ iVignali et a l.' '2010^, Spitzer 
(Zakamska et al. 2008), spectropolarimetry (Zakamska et alJ 
2005) Gemini ( Liu et al. 2009), and the VLA dZakamska et al] 
,2004c iLal & Ha201 0) yield a broad view of the properties of 
the optically selected obscured quasar population. We target the 
low-redshift end of the sample, to maximize our spatial resolu- 
tion of the NLR. I n our first paper, we examined the host galax- 
ies of our targets (iGreene et al. I I2OO9I Paper 1 hereafter). Here 
we study the spatial distribution and kinematics of the ionized 
gas. After describing the sample and observations (§2), we turn 
to the NLR sizes (§3) and then the spatially resolved kinematics 
of the sample as a whole (§4). We present two candidate dual 
obscured AGNs (§5) and then summarize and conclude (§6). 

2. THE SAMPLE, OBSERVATIONS, AND DATA REDUCTION 

The sample and data reduction were introduced in detail in 
Paper 1 (Table 1). The sample was selected from iReves et al.l 
(2008). We focused on targets with z < 0.45 to ensure that 
[O III] A5007 was accessible in the observing window and im- 
posed a luminosity cut on the [O III] line of L[om] > lO"'^ erg s"' 
to pre-select luminous quasars (estimated intrinsic luminosity 
Mb < -24 mag). Radio flux densities at 1.4 GHz were ob- 
tained from the Faint Images of th e Radio Sky at Twenty cm 
survey (FIRST; Becker et al."1995') and the NRAO VLA Sky 
Survey (NVSS; .Condon et al...l998.) . With one exception, aU 
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Spatial Position (arcsec) Spatial Position (arcsec) Spatial Position (arcsec) 

Fig. 1. — Three example [O iii] (solid) and continuum (dotted) surface brightness profiles as a function of spatial position, plotted in units of 
10~" erg s"' cm"^ arcsec"^. In each case the line spectrum was summed over a velocity twice that of the FWHM of [O iii], while the continuum was summed over 
a Une-free band 10 times wider The inferred [O iii] (short-dashed line) diameters ai'e shown for reference. In the center panel, we include a flux calibration star for 
reference (dashed). 



objects are radio-quiet, as dete rmined by their position on the 
Lronii -i^L^ (1.4 GHz) diagram (!Xu et al.ll999l:IZakamska eTal] 
|2004|) . and they are at least an order of magnitude below the 
nominal radio-loud vs. radio-quiet separation line in this plane. 
The single radio-loud object in the sample, SDSS J1124H-0456, 
is a double-lobed radio galaxy (alternate name 4CH-05.50) with 
vLi, (1.4 GHz)= 2.9 X 10"^^ erg s~' which was observed with 
a slit nearly perpendicular to the orientation of its large-scale 
radio lobes. 

We observed 15 objects over two observing runs us- 
ing the Low-Dispersion Survey Spectrograph (LDSS3; 
lAUington-Smith et al.lll994 ') with a 1" x 4' sHt at the Magel- 
lan/Clay telescope on Las Campanas. The seeing was typically 
^1" over the two runs. We integrated for at least one hour per 
target and covered one or two slit positions (Table 1). Lower-z 
targets were observed with the VPH-Blue grism in the reddest 
setting, for a wavelength coverage of 4300 - 7050 A, while the 
higher-z targets were observed with the bluest setting of the 
VPH-Red giism (5800-9400 A). The velocity resolution in 
each setting is ainst ~ 67 km s"'. In addition to the primary 
science targets, at least two flux calibrator stars were observed 
per night and a library of velocity template stars consisting of 
F-M giants was observed over the course of the run. 

Since we have only long-slit observations, we do not sample 
the full velocity field of the gas or stars in the galaxy. With 
a few exceptions, the galaxy images were only marginally re- 
solved in the SDSS images. Thus in selecting position an- 
gles to observe we were mainly guided by visual inspection 
of the color composite images. Since these galaxies typi- 
cally have very high equivalent width [O III] lines, we at- 
tempted to identify [O III] structures based on color-gradients 
in the images. As a result, the slit is not necessarily oriented 
along the major or minor axis of a given galaxy. In partic- 
ular, it is important to keep in mind when judging the ra- 
dial velocity curves of the spiral galaxies (SDSS J1106H-0357, 
SDSS J1222-0007, SDSS J1253-0341, SDSS J2126H-0035 and 
likely SDSS J1124H-0456). Of these, SDSS J1106H-0357 and 
SDSS J2126H-0035 were observed along the major axis, and 
SDSS J1222-0007 is within 20° of the major axis. The others 
are observed at ^ 40° from the major axis. None were observed 
solely along the minor axis. 

Cosmic-ray removal w as performed usin g the spectroscopic 
version of LACosmic (Ivan Dokkumll200l . and bias subtrac- 



tion, flat-field correction, wavelength calibration, pattern-noise 
removal (see Paper 1), and rectification were performed using 
the Carnegie Observatories reduction package COSMOSQ. For 
the two-dimensional analysis discussed in this paper (e.g., the 
[O III] size determinations) we additionally use the sky subtrac- 
tion provided by COSMOS. The flux calibration correction is 
determined from the extracted stan dard star u sing IDL routines 
following methods described in Mat heson et al. (2008) and then 
applied in two dimensions. In the first paper we demonstrate 
that the absolute normalization of the flux calibration is reliable 
at the ^ 40% level. "Nuclear" measurements refer to the 2. "25 
spatial extraction. 

3. NEBULAR SIZES 

3.1. Measurements 

The physical extent of the NLR provides one basic probe of 
the impact of the AGN on the suiTounding galaxy. We work 
with the rectified two-dimensional spectra. In order to boost 
the signal in the spatial direction, we collapse each spectrum 
in the velocity direction. We use a band with a velocity width 
that is twice the full width at half-maximum (FWHM) of the 
nuclear [O III] and centered on the nuclear [O III] line (Fig. 
[TJ. The line width is measured from a continuum-subtracted 
spectrum, but we do not perform continuum subtraction on the 
two-dimensional spectra. This high signal-to-noise (S/N) spa- 
tial cut allows us to measure the NLR sizes much more sensi- 
tively than from typical narrow-band imaging. Specifically, we 
measure the total spatial extent of the line emission down to a 
5 E limit, where E is determined from spatially-offset regions 
of the collapsed surface brightness profile. We are reaching typ- 
ical depths of E « 10"'^ erg s"' cm~^ arcsec"^. In three cases 
the nebular spectra are not spatially resolved (i.e., the spatial 
distribution matches that of a standard star). There are six ob- 
jects for which we have multiple slit positions. The range in 
nebular size derived from cases with multiple slit positions is 
~ 30%. 

In a few cases (SDSS J1356-1026, SDSS J2126+0035 & 
SDSS J2212-0944), the line ratios change as a function of ra- 
dius and [O III]/H/3 falls below three. This changing ratio may 
reflect changes in the ionization parameter or gas-phase metal- 
licity, or a transition from ion ization dominated by the AGN 
to H II regions (e.g., iBennert et al. 2006). By ionization pa- 
rameter, we mean the ratio of the density of ionizing photons 
to the density of electrons. Given the luminosities of quasars 
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40 42 44 3 3.5 4 2 2.5 [ 

log (L[o ,„] / erg s-i) log / pc) log (FWRMf^ „„ / km s"') 

Fig. 2. — Radii of the [O iii]-emitting regions, as measured from spatial cuts through the spectra centered on the [O iii] line, compared with various other properties. 
Our measurements ai'e shown in large circles and upper limits on the nebul ar sizes, in u nresolved cases, are shown as aiTo ws. Our estim ated uncertainty of ~ 3 is 
demonstrated f or one point only. We also show measureme nts from .Bennert et al.H2002l broad-hne AGNs; filled ?riang/e^). lSchmitt et aU<2Q03. grey filled squares), 
IFraguelli et al.l f2003. asterisks), and lBennert et al.l j2006l crosses), (a): i?[oiii) compared to total L[oin]- The two quantities are correlated (Kendall's r = 0.89 with 
probability P < 10"^ that no correlation is present). Our best-fit is shown with a solid line. Our emission-line regions are considerably larger at a given L[oiiij than 
either the broad-line quasars considered by Bennert et al. (2002; filled triangles and dotted line) or the relation derived by Schmitt et al. (2003; dashed line), likely 
due to differences in depth. Thus, we exclude the latter two in our fitting, (b): R[om] compared to the size of the galaxy continuum measured in a comparable fashion 
(see text). Here we include only the sample of Bennert et al. (2006) for comparison, since we do not have galaxy sizes available in the other cases. In general it is 
clear that the continuum and emission line sizes are of the same order (Kendall's t = 0.52 with P = 0.098). The solid line shows the 1 : 1 relation to guide the eye. 
(c): S[oiiij plotted against the FWHM of the [O iii] line from the nuclear extraction. Here we include FWHM measurements from Whittli^ il992al) for the Bennert et 
al. (2006) and Schmitt et al. (2003) samples. There is clearly a trend (Ken dall's r = 0.75,P = 3 X 10^). The two oufliers (bottom, right) are Mrk 3 and NGC 1068. 
However, we note that the size of NGC 1068 is likely underestimated here; i Veilleux et al] J2003I) find a size of ~ 11 kpc using a tunable filter 



in our sample and the rates of star formation in their hosts 
jZakamska et alj l2008l) . we expect that the number of ioniz- 
ing photons from the quasars exceeds that from stars by about 
an order of magnitude. Nevertheless, since quasar illumination 
is not necessarily isotropic and since photons from star forma- 
tion are distributed more uniformly within the galaxy than those 
arising from the central engine, it is plausible that we may see 
gas excited by stars in the outer regions of the galaxy. Shock 
excitation is unlikely since the linewidths are uniformly narrow 
in these outer regions. To be safe, we exclude the regions with 
[O III]/H/3< 3 when calculating the NLR sizes. We have not 
applied any correction for redde ning, which c ould be substan- 
tial (e.g., Melendez et al. 2008). iReves et al] (2008) show that 
deriving robust extinction corrections for the SDSS obscured 
quasars is not straightforward, and we neglect such corrections 
here. 

One of the objects in our sample, SDSS J1356-I-1026, has 
a much more dramatic extended emission-line nebula than the 
rest (Figure|3]l. We will discuss the detailed kinematics and en- 
ergetics of this object in more detail in a parallel paper (Greene 
& Zakamska in preparation). For the present work, we explore 
the implications of detecting one single extended emission-line 
region in the sample (§6). 

We should note that deriving nebular sizes is an ill-defined 
task. First of all, ionized nebulae need not have regular shapes, 
and so the definition of size is not necessarily well-defined. 
This difficulty is only exacerbated when long-slit spectra are 
used to define the size, since our slit may well miss spatially 
extended regions. Furthermore, the concept of size depends 
sensitively on the depth of the observation. Deep observations 
probing depths of a few times 10""^ erg s"' cm"^ arcsec"^ in- 



deed reveal faint, extend ed gas with a range of morphologies 
(e.g., IVeilleux et al.ll2003[) . Thus, the primary size uncertainties 
are in these systematics, which dwarf the measurement errors. 

We quantify the uncertainties in the following manner First 
of all, since we have not included surface brightness dimming, 
there is a dispersion of ~ 0.1 dex in the sizes due to dis- 
tance. Secondly, and more important, narrow-line regions are 
not strictly round. Thus, depending on the position of the long- 
slit, we may derive a different answer We have found, for 
the six objects with multiple slit positions, that the sizes agree 
within 30%. Finally, and most difficult to quantify, the shape 
will likely grow more irregular as we push to lower flux limits . 
We have attempted to quantify this dispersion in emission line 
profile using the ratio between the luminosity-weighted mean 
width of each spatial [O III] profile and the adopted radius mea- 
sured down to a fixed surface brightness. Were the NLRs all of 
the same shape, then the mean width would be a fixed fraction 
of the total size. Instead, the ratio ranges from 0.2 to 4, with 
a typical value of three. We thus adopt a factor of three as the 
overall uncertainty in the sizes. Additionally, we flag as particu- 
larly uncertain those systems with a nearby massive companion 
galaxy, since there we are further contaminated by tidal gas. 

3.2. Obscuration, Ionization, and Excitation 

N LR sizes have been measured from narrow-band imag- 
ing (" Mulchaev et al.lll996al: iBennert et an i2002'; Schmitt et d] 
2003 ) and from long-slit spectroscopy (e .g., lUngeret al.J987l : 
Fraq uelli et al.ll2003l; iBennert et al.ll2006h . Narrow-band imag- 
ing is preferable for studying the NLR morphology, but reaches 
shallower limits than the spectroscopy. 
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Fig. 3. — Left: Image (r-band) of the merging galaxies SDSS J1356+1026. 
Scale bar indicates 10" (22 kpc). N is down and E to the right. Right: The 
two-dimensional spectrum centered on the continuum-subtracted [O iii] A5007 
line. Oriented as the image. The spectrum spans 1814 km s^'in the velocity 
(x) direction. 

Integral-fiel d observations allow on e to study two-dimensional 
kinematics (iHumphrev et a l. 2010), but for local objects only 
cover the inner NLR (e.g. jBarbosa et al- 2009). 

We compile a comparison sample of lower-luminosity ob- 
scured AGNs with meas ured NLR size s from the literature 
dBennert et alj 12002; Sch mitt et all 120031: iBennert et all l2006t 
iFraquelli et ailHoOSI) . We include the Bennert et al. (2002) 
and Schmitt et al. (2003) measurements in Fig. 2 for com- 
pleteness, but note that the sizes cannot be compared directly 
with those we measure here, because of the difference in depth. 
The limiting surface brightness values that we achieve in this 
work are at least a factor of 10 deeper than these narrow -band 
imaging studies from space, which range from ~ 10"'^ -3 x 
10"'"^ erg s"' cm"^ arcsec"^. For this reason, we do not include 
the space-based measurements in any analysis presented here 
(e.g., fitting of relationships). Fraquelli et al. do not quote 
sizes but rather provide power-law fits to the surface bright- 
ness as a function of distance to the nucleus. Taking their func- 
tional form, we calculate sizes that match our limiting surface 
brightness of ^ 10"'^ erg s"' cm"^ arcsec"^. For uniformity, 
we calculate sizes for Bennert et al. (2006) in the same way, 
and we adopt their smaller radii in cases where star forma- 
tion dominates in the out er parts. I n cases of overlap between 
works, we prefer the Ben nert et al.l (2006) observations, since 
they are both sensitive and take into account photoionization by 
starlight. The measurements for our sample are summarized in 
Figure|2]and Table 2, while the comparison samples are shown 
in Figure |2] 

The observed distribution of NLR gas depends on the geom- 
etry and luminosity of the ionizing source, the geometry and 
kinematics of the host ISM (e.g., disk, spherical, outflow, or 
infall), and the density distribution of the gas. While most of 
these are likely related to the morphology and dynamical state 
of the galaxy, the geometry of the ionizing source is tied to the 
orientation of the AGN. In the simplest model, the galaxy ISM 
is spherically distributed, while the ionizing radiation from the 
AGN emerges anisotropically along lines of sight unaffected 
by the circumnuclear 'torus', as postulated by unified mod- 
els of AGN activity. In this case we expect to see ionization 
cones when the beam is not pointed directly at us, reflecting 
the geometry of circumnuclear obscuration. Su ch cones are 
observed in images of nearby Seyfert galaxies (|Poggel[T988l: 



iTadhunter & Tsvetanovl [T989t IStorchi-Bergmann et al.l [19921) 
and more recently in the l uminous obscured quasars studied 
here (IZakamska et al.ll2006h . 

In this simplest geometry, we would expect to find smaller 
sizes in unobscured sources, when looking closer to the axis of 
the ionization cones. The difference in distributions depends on 
the expected opening angle of the torus, with larger difference 
for smaller opening angles. Recent observations of obscured 
quasars suggest that the sp ace densities of ob scured and unob- 
scured sources are ~ equal dReves et al.l2008l) . leading to open- 
ing angles of < 120°, but even if significantly smaller opening 
angles are assumed, the expected differences in the median pro- 
jected size between the two populations is small (< 0.2 dex). At 
low redshift and (thus) lower luminosity, io nization cones are 
also observed in unobs cured sources (e.g., lEvans et al1ll993l : 
iBokse nberg et al.' '1995'). while round NLRs are observed in 
both types (Mulchaey et al. 1996b). Presuma bly, the ISM is 
not always spherically dist ributed or relaxed (iMulchaev et al.l 
II996bHSchmidt et alJl2007 l). 

In Figure I2J2 it appears that the NLR sizes of our obscured 
quasars are larger than the unobscured ones (median difference 
0.4 dex). However, this difference can be explained by differ- 
ences in the depths of the observations, thus we cannot address 
orientation differences in detail from these samp les. It i s inter- 
esting to note that at lower luminosities, Schmit t et al. I (I2003h 
do not see a significant size difference between the two popula- 
tions. These observations, while shallow, are uniform between 
the obscured and unobscured populations. 

There has been some debate in the literature about the 
slope of a purported corre lation between the NLR size and 
the AGN luminosi ty (e.g ., iBennert et all 120021; ISchmitt et all 
120031; IBennert et al.ll2006l) . Some correlation is expected, given 
that the AGN is photoionizing the NLR gas, but the form it 
takes may tell us something about covering factor or density 
as a function of luminosity. It is clear from Figure |2jfl) that 
generally larger NLRs are found in more luminous objects 
(Kendall's r = 0.89 with probability P < 10"^ that no corre- 
lation is present). It is also clear that there is substantial scat- 
ter; we find an rms scatter of 0.3 dex in radius at fixed L[oiii]. 
We performed Monte Carlo simulations of ionization cones ob- 
served at random directions (restricted to be outside the cones). 
They suggest that the orientation of the NLR axis relative to the 
line of sight is not a significant source of the observed scatter. 
At a fixed NLR size, orientation effects introduce a scatter of 
< 0.15 dex within each (obscured or unobscured) subpopula- 
tion, even when a wide range of opening angles is allowed for. 
Therefore, the observed scatter is likely due to the combination 
of the true variance in NLR sizes at a given luminosity and to 
the difference s in the definition of NLR "size". For instance, 
IBennert et al.l (12006.) derive sizes that are factors of ^ 2 larger 
than those based on HST narrow-band imaging because of their 
increased sensitivity. Given that the NLR is not always spher- 
ically symmetric or smooth, defining a meaningful size that is 
insensitive to depth is a difficult problem. 

For completeness, we fit a power-law relation between L[oiii] 
and NLR size, using all narrow-line comparison samples as 
well as the objects considered here. Because there are up- 
per limits on the sizes, we calculate a linear regression us- 
ing the binned Schmitt method, from the Astronomy Sur- 
vival Analysis (AS URV) software as imple mented in iraf 
dFeigelson & NelsonlllHlllsobe et alJll986l) . The fit is shown 
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Fig. 4. — Resolved spectroscopic measurements for SDSS J1222-0007. In the top panel we show the central spectrum, extracted within the 
inner ~ l", arbitrarily normalized. Then we show the flux in the [O m] line (filled black squares) and the continuum (open stars) as a function of 
radius (10"'^ erg cm"^ s"'). In the third panel we show the radial velocity curve of the peak of [O m] (open triangles), the velocity of the peak of 
the H/3 line (open circles), and the luminosity-weighted mean velocity of [O m] (open squares). In the fourth panel the velocity dispersion of the 
[O III] line (open black circles) is compared with the maximum velocity (the velocity at 20% of the line maximum) to the red (red open circle with 
cross) and the blue (blue open circle with dot). In the final panel the ratio of [O iii]/H/3 is plotted (black crosses). A ratio of [O in]/H/? of three is 
noted with the dashed line. 
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in Figure 2. We find: 

log(/;NLR/10V) = (0.22 ± 0.04) log(L[o„i]/1042gj.g ^-i) 

+(3.76 ±0.07). 

The shallow slope we observe is consistent with a picture in 
which the nebulae are matter-bounded. At the distances from 
the quasar that we are probing with our observations, the den- 
sity of material is low enough that the emissivity is no longer 
limited by the flux of photons by the quasar, but rather by the 
low density of the gas, and a large fraction of photons can es- 
cape into the intergalactic medium. Note that the correlation 
between AG N continuu m luminosity and L[oiu] in broad-line 
AGNs (e.g.. [19801) suggests that the nebulae are limited 
by the number of photons in the bright central regions of the 
galaxy, but that the situation changes in the diffuse outer parts. 
If so, we would expect size to scale as the square-root of lumi- 
nosity at low luminosities and then flatten out to at high lumi- 
nosities, modulo differences in host galaxies. 

In addition to measuring the nebular sizes, we also parame- 
terize the luminosity drop in the outer parts as a power-law and 
measure the power law slope (L[oiu]Oc r~^). The slopes range 
from 1 .6 < /? < 3.8 (Table 2). These slopes correspond to den- 
sity profiles with slopes ranging from 1.3 to 2.4, in good agree- 
ment with the HST observations of Zakamska et al. (2006). 

One concern, as pointed out by Netzer et al. (2004), is that 
eventually the ionizing photons will run out of interstellar 
medium to ionize, particularly in the most luminous quasars. 
The NLR size cannot in general grow indefinitely beyond the 
confines of the host galaxies. In Figure |2/7), we compare the 
continuum and nebular sizes. Rather than using effective radii 
of host galaxies from photometry, we use the same method to 
measure the continuum extent as we used for the [O III] lines, 
collapsing the two-dimensional spectrum in the spectral direc- 
tion over line-free regions to boost the signal. Galaxy sizes 
are weakly correlated with NLR size (Kendall's r = 0.52 with 
P = 0.0975). We see that the NLR sizes are comparable to the 
galaxy continua. The exception is SDSS J1356+1026, which 
contains the spectacular bubble shown below. At these lumi- 
nosities, the AGN is effectively capable of photoionizing the 
entire galaxy ISM, as well as companion galaxies out to sev- 
eral tens of k pc, as we saw in some of our previous long-slit 
observations (ILiu et al.ll2009l) . 

Outflowing com ponents of the NLR are routinely seen in ra - 
dio galaxies (e.g., (McCarflivl ll9 93t IVillar-Martfn et al.lll999h 
as weU as in Seyfert galaxies (e.g. .^ Crenshaw & Kraemed2000t 
iRupke et al. l 2005). On small scales, detailed modeling of the 
inner (< 300 pc) NLRs of a few local AGNs with HST indi- 
cates a surprising uniformi ty in behavior, with v c>c r along an 
evacuated bicone (Crensha w & Kraemeill2000HCrenshaw et al.l 
l2000t iRuiz e"tani2001 ). Interestingly, we see similar qualita- 
tive behavior in SDSS J1356+1026 (Greene & Zakamska in 
prep). However, in general, NLR kinematics on larger scales 
are not as uniform, with mechanisms rang ing from jet ac- 
celeration to radiation pressure drivin g (e.g., iRuiz et all 1200 it 
iGroves et al."2004'; 'RuokeetaDHop). For a complete review, 
see Veilleux et al. ( 2005). We do find some correlation between 
FWHM and luminosity (Figure |2}:; Kendall's r = 0.75, P = 
3 X 10"^). We will argue below based on the observed large 
velocity dispersions at large radius that the AGN energy is stir- 
ring up the gas on large scales, thus explaining this correlation. 

4. TWO-DIMENSIONAL ANALYSIS 
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Fig. 5. — Maximum velocity as a function of radius. The maximum velocity 
is measured as the velocity at 20% of the line peak, as referenced to the sys- 
temic velocity of the stars. At each position we show the maximum red- (red 
circle) and blue-shifted (blue circle) velocity. While there is extended emission 
in nearly all targets, the typical velocities are not higher than ~ 400 km s"' for 
the majority of the targets. SDSS 11356+1026 is highlighted with large black 
circles. The dashed line highlights a velocity of 500 km s"', which is the ap- 
proximate escape velocity for these galaxies at r^. We note that the effective 
radii of these galaxies, for which we have well-resolved imaging, range from 
1 lb 0.4 to 1 1 lb 2 kpc, with a median of 3 ib 0.6 kpc. 

In this section we present the results of our two-dimensional 
analysis on the long-slit spectra. First we present velocity and 
dispersion profiles, as well as emission line ratios, as a function 
of position. To obtain these measurements, we extract spectra 
at uniform intervals as a function of spatial position along the 
slit. We start with rectified two-dimensional spectra from COS- 
MOS. Each spectrum is extracted with a width of 0"95 (5 pix- 
els) to match the typical seeing of the observations. The central 
spatial position is determined by the spatial peak in the [O III] 
emission. The systemic velocity is determined from the absorp- 
tion lines. Galaxy continuum subtraction is performed for each 
spectrum using a scaled version of our best-fit model from the 
nuclear spectrum, with only the overall amplitude allowed to 
vary. While this is not strictly speaking a correct model, we 
have insufficient S/N in the off-nuclear spectra to constrain ve- 
locity or velocity dispersion, let alone changes in stellar popu- 
lations. 

Once the continuum-subtracted spectra are in hand, we fit 
the H/3h-[0 III] AA4959, 5007 lines for each spectru m as in Pa- 
per I (see also iHo et al.l[T997l iGreene & Holl2005l) . Each line 
is modeled as a sum of Gaussians (a maximum of two for H/3 
and three for [O III]). The relative wavelengths of each transi- 
tion and the ratio of the [O III] lines are fixed to their laboratory 
values, but the central velocity and line widths are allowed to 
vary from spectrum to spectrum. From these fits we are able to 
derive velocity, velocity dispersion, and line-ratio profiles as a 
function of spatial position. We report three measures of veloc- 
ity at a given position, the peak in the [O III] line, the peak in 
the H/3 line, and the flux-weighted mean velocity in the [O III] 
line. The velocity dispersion is measured as the FWHM of the 
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[O III] model divided by 2.35. At each spatial position we also 
measure the "maximum" and "minimum" velocitie s as the ve- 
lociti es at 20% of the [O III] peak intensity (e.g., Rup ke etal] 
I2OO2I) relative to the systemic velocity of the stars (shown as 
blue bullseyes and red crosses, respectively, in Figure 0) 

Errors are derived from Monte Carlo simulations. For each 
spectrum we generate 100 mock spectra using the best- fit pa- 
rameters at that radius and the S/N of the original spectra. We 
fit each mock spectrum and the quoted parameter errors encom- 
pass 68% of the mock fit values. 



4. 1 . Radial Velocity Curves 

In Figure |4] we present a representative radial velocity 
curve for SDSS J1222-0007. The remainder are shown in 
the Appendix. First, we note that overall the radial ve- 
locity curves are flat. In Paper I we presented detailed 
two-dimensional photometric fitting of these galaxies (with 
the exception of SDSS J1124H-0456 and SDSS Jl 142H-1027). 
Using these fits, we divide the sample by the bulge-to- 
total ratio (B/T), and call galaxies with B/T< 0.1 disks 
(SDSS J1106H-0357, SDSS J1222-0007, SDSS J1253-0341, 
and probably SDSS J1124H-0456), while the rest are bulge 
dominated. Additionally, those with clear tidal signa- 
tures are "disturbed" (SDSS J0841H-0101, SDSS J1222-0007, 
SDSS J1356H-1026, and SDSS J2212-0944). We would expect 
to see the signature of rotation most clearly in disk-dominated 
galaxies. We note once again that SDSS J1106H-0357 and 
SDSS J2126H-0035 were observed along the major axis, 
SDSS J1222-0007 was within 20° of the major axis, and the 
remaining two galaxies were observed at a ^ 40° angle to the 
major axis. We would expect to see the signature of rotation in 
most of these galaxies. Instead, we only see rotation in the case 
of SDSS J1106H-0357, SDSS J1124H-0456, SDSS J1142H-1027, 
and SDSS J2212-0007. Although with such a wide range of 
position angles, and such a small sample, it is hard to say for 
sure, we find it suggestive that neither SDSS J1253-0341 nor 
SDSS J2126H-0035 shows rotation. 

The sample galaxies showing rotation in their radial veloci- 
ties also tend to show declines in Cg by factors of two or more 
in the outer parts (e.g., SDSS J1124H-0456). In contrast, those 
galaxies with flat radial velocity curves (the majority in this 
sample) also have notably flat dg distributions at kpc scales. 
Again, this is strongly in contrast to the ki nematics in the 
stars, even in bulge-dominated systems (e.g., iJorgensen et al.l 
Il995h . More to the point, it is in contrast to the kinemat- 
ics of warm gas in inactive late-type (e.g., Pizzella et al. 20041) 
and early-type (e.g., Fillm ore et al. 1986; Bertola et al. 1995t 
IVega Beltran et al. 12001 1) spiral galaxies. 

In Paper I we showed that Cg in the nucleus is uncorrected 
with (T* . Again, this behavior is in striking contrast not only 
to inactive galaxies but also to local, lower-luminosity active 
galaxies, for w hich it has been long been known that on average 
CTg/2*« 1 (e.g.. lNelson & Whittle! 1996tlGreene & Hol2005tlHol 
Here we are making a stronger statement. Not only is 
the luminosity-weighted gas dispersion uncorrelated with the 
dispersion in the stars, but the dispersion in the gas stays high 
out to kpc scales in these galaxies. These observations provide 
new reason to doubt that gas velocity dispersions can be sub- 
stituted for s tellar velocity dispersions in luminous AGNs (e.g., 
[Shields et al.ii2003;.Salviander et al...2007.) . 
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Fig. 6. — Nucleai' [O iii] luminosity compai'ed witli tiie maximum gas veloc- 
ity for each object. The maximum velocity plotted here is the largest that we 
measure for any individual object, although we exclude the maximum velocity 
in the nuclear spectrum. No correlation is apparent between the two (Kendall's 
r = 0.13 with a probability, P = 0.73 of no correlation). 

This behavior is different from that seen in regular inactive 
galaxies. It is also different from that in local, well-observed 
Seyfert galaxies. Previous work looking at the kinematics of 
lower-luminosity local Seyfert gal axies has found ev idence for 
a two-tiered NLR structure (e.g., lUnger et"arill987l) . In such 
objects, the inner or classical NLR extends to a few hundred 
pc and has linewidths of > 500 km s~' . At higher spatial res- 
olution, there is clear evidence for out flow in the inner hun- 
dreds of pc in well-studied objects (e.g..lCren shaw et al."2 000l; 
ICrenshaw & Kraemed l2000t iRuiz et al. iMoi; Barbosa e t al.l 
'2009^. In contrast, at larger radius, the linewidths drop and 
the kinematics of the NLR gas simpl y reflect that of the bulge 



or disk in which the gas sits (e.g., ISfelson & Whittle! 119961; 



iGreene & H0II2OOI IWalshetal.1 120081). Clearly, the observed 
kinematics of gas in hosts of obscured quasars are quite dissim- 
ilar from this picture. 

Many of the host galaxies of our obscured quasars have 
nearby companions and/or show signs of recent interactions. It 
is therefore possible that the gas is being stirred by gravitational 
interactions with nearby galaxies. To explore that possibil- 
ity further, we examine the analogous inactive ul tra-luminous 
infrared galaxies. The integral-field spectra of IColina et al.l 
(j2005) show that even in these ongoing mergers the gas kine- 
matics traces that of the stars. The cTg profile is typically seen 
to decline in the outer parts as in non-merging systems, again 
in contrast to our findings for hosts of obscured quasars. Of 
course, there are exceptions in the Colina et al. sample, where 
the gas velocity dispersions are very complex. On the other 
hand, the mergers are more advanced in general than in our 
sample. Thus, while we cannot rule out gravitational effects in 
all cases, it seems most likely that the nuclear activity is directly 
responsible for stirring up the gas. We now address whether 
there is evidence for bulk motions (e.g., large-scale outflows) 
in the gas based on the kinematics. 
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Fig. 7. — Fraction of the [Oiii] emission line luminosity emerging at a veloc- 
ity greater than the escape velocity for all objects where the escaping fraction 
exceeds 10% in at least one shell. Object names are labeled with the first 
four digits of the full position. The escape velocity is approximated from the 
measured velocity dispersion assuming that vv- = \/2i7» and that the halo has a 
radius of 100 kpc. We then simply integrate the part of the line with velocities 
greater than the escape velocity. Here we show the fraction of the line emis- 
sion that is in this high- velocity component as a function of radius from the slit 
center. 



4.2. Maximum Velocities 

We have derived "maximum" red- and blueshifted veloci- 
ties at 20% of the line profile, relative to the systemic veloc- 
ity of the stars. We examine the distribution of maximum 
velocities as a function of radius for the ensemble of spec- 
tra in Figure |5] While the emission extends to kpc scales 
for the majority of the targets, the gas velocities are not typ- 
ically very high. The median maximum blue velocity at 8 
kpc is (vbiue) = 400 ± 70 km s"' while towards the red it is 
(i^'red) = 330 ±50 km s"', where we quote errors in the mean. 
A few objects (SDSS J1253-0341, SDSS J1222-0007) have 
gas at velocities exceeding 500 km s"'. We note that the ef- 
fective radii of these galaxies, for which we have well-resolved 
imaging, range from 1 ± 0.4 to 11 ± 2 kpc, with a median of 
3 ± 0.6 kpc. These velocities exceed the velocity dispersions 
of the galaxies, but they do not compare to the ^ thousands 
of km s"' outflow velocities seen by Tremonti et al. (2007) and 
postulated to be driven by recent AGN activity. Furthermore, 
they are not close to the escape velocity needed to actually un- 
bind the gas. As we show in Figure|6] there is no evidence for a 
correlation between the nuclear L[om] luminosity and the maxi- 
mum observed velocity (Kendall's r = -0.10 with a probability 
f = 0.8 of no correlation). 

We now quantitatively address whether any of the gas is ap- 
proaching the escape velocity. Following Rupke et al. (2002), 
we calculate an approximate escape velocity for each galaxy by 
assuming that the circular velocity scales with the velocity dis- 
persion as vv = -\/2cr* . Assuming the potential of an isothermal 
sphere, the escape velocity as a function of radius scales as: 

Vesc(r) = \/2y, [ 1 + ln( 1 + w/r)]"'' . (2) 

Although Tmax IS unknown, the escape velocity depends only 
weakly on its value. Thus we assume rmax =100 kpc in all 



cases. The escape velocities thus estimated range from 500 to 
1000 km s"' over the entire sample, but only vary by ^ 60% for 
an individual object over the range of radii that we probe. 

With escape velocities in hand, we can now address what 
fraction of the line emission comes from gas that is moving at 
or above the escape velocity. We first ask whether there is gas 
exceeding the escape velocity at each radius. With the same 
definition of systemic velocity as above, we integrate the line 
emission that exceeds the escape velocity to either the red or 
blue side of the systemic velocity. We then normalize by the 
total flux at that radius. These fractions are plotted as a func- 
tion of radius in Figure |7]for the 5 objects in which at least 10% 
of the gas is nominally escaping for at least one radial position. 
For illustrative purposes, we focus here on the blue-shifted gas. 
In addition to calculating the escaping fraction at a given radius, 
we can also calculate an overall escaping fraction. They range 
from < 1 % - 25% with a median value of 2%. 

Nominally only a small fraction of the NLR gas is moving 
out of the galaxy at or around the escape velocity. However, the 
projection effects may be severe, and especially so because in 
obscured objects the gas motions are expected to occur largely 
in the plane of the sky. Therefore, our estimates are a lower 
limit on the actual escape fractions (see §6 for details). Further- 
more, as discussed further below, we have good reason to think 
that the medium is clumpy. Depending on whether the out- 
flowing component has the same clumping factor as the bound 
gas, it is difficult to translate these observed fractions into mass 
fractions. 

In addition to the escaping fraction, we would like to know 
how much mass is involved in the outflow. The standard 
method of estimating the density of the emission line gas uses 
density diagnostics such as the ratio of [S II] A6716/A6731 
or [O II]A3729/A3726. Neither of these is available in the 
Magellan spectra, and with several hundred km s"' velocities, 
the [O II] doublet is blended enough to be difficult to mea- 
sure. The continuum-subtracted SDSS spectra that integrate 
all emission within the 3" fiber yield a measurement of the 
[S II] A6716/A6731 ratio for all but the highest redshifts. Us- 
ing the iraf task temden, these can be translated into den- 
sities ranging from 250-500 cm"-', with a mean of 335 cm"-'. 
These values are consistent with those commonly seen in spa- 
tially resolved observat ions of extended NLRs and used in 
mass estimations (e.g., iNesvadba et al.l 120061; iFu & StocktonI 
2009 and many others). However, such measurements can be 
highly biased toward high densities in clumpy gas. Specifi- 
cally, the recombination line luminosity depends on density as 
L cx J dVarietip, whereas mass goes like M cx J dVnp, so the 
mass of the gas, its density and degree of dumpiness and its 
line luminosity are related through 

Mg = 1 .7 X lO'^L^'^ («,)"' k"'M0 . (3) 

Here we used a recombination coefficient a = 1.62 x 
10"'"^ cm^ s"' appropriate for a 20,000 K gas, and k = (n^) /{rie)^ 
is the degree of dumpiness, which by definition is > 1 and 
can be substantially greater We have adopted a higher tem- 
perature than typical based on the [O III] A4363/[0 III] A5007 
line ratio. In general, the ratio ranges from ^ 0.01 in the cen- 
tral regions to ^ 0.04 further out. These ratios correspond to 
T = 11000-23000K, and thus we adopt a temperature repre- 
sentative of the outer regions. 

The standard method of calculating the mass involved 
amounts to using this equation with k = 1 and n^. of a few x 100 
cm"-', and produces an absolute minimum on the gas mass vis- 
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Fig. 8. — The measurements of the scattered light emission made in the rest-frame UV band using the HST can be translated into a constraint on the density of 
scattering particles. Grey points are all measurements from the HST images of Zakamska et al. (2006), while fitting was done with the large symbols, and open 
symbols are binned. The observed surface brightness of scattered light (vertical axis) as a function of distance r from the center of the galaxy (horizontal axis) is 
proportional to the illuminating flux from the central AGN, L/iAnr^), to the scattering cross-section of particles, dcr/dfi, and to the column density of scatterers 
~ n,r/3, where /3 is the opening angle of the scattering region. Solid lines show the best power-law fit to electron density necessary to reproduce the scattered 
emission, assuming that scattering is due to dust particles. Per given hydrogen density, dust particles are about 60 times more efficient scatterers than electrons, even 
if the gas is fully ionized. The measurements are made using host-galaxy-subtracted images at about rest-frame 3100A. The deficit of light in the central parts is due 
to the PSF smearing and dust obscuration. 



ible in the emission lines of a few xIO^Mq. However, such 
high densities are in direct conflict with our observations. For 
one thing, we see high [O III]/H/3 ratios, and thus high ion- 
ization parameters, and presumably low densities, at large ra- 
dius. Also, the observed extended scattering regions in ob- 
scured quasars place an independent constraint on gas densi- 
ties dZakamska et aL 2006). Scattered light flux is oc J dVn.,, 
where is the density of scattering particles, electrons or dust. 
Assuming purely electron scattering, HST observations can be 
fit by density profiles that decline as ^-(i 5-2.5) \f^ith den- 
sity (wp) = 1 cm"-' at a distance of about 3 kpc from the cen- 
ter (Figure [Sj. The scattering angle is not well known, but 
it introduces only about a factor of two uncertainty in this 
measurement. Dust particles are even more efficient scatter- 
ers than electrons, so in the more realistic case of dust scatter- 
ing, w hich is suggested by several lines of observational evi- 
dence JZakamska et al.ll2005i) . the implied mean density is con- 
strained to be even smaller, («(,)(lkpc) ^ 0.016 cm"^. The un- 
certainties are larger in the case of dust scattering, because the 
density measurement is sensitive to the assumed gas-to-dust ra- 
tio and the scattering angle ( for this parti cular value, 90° and 
Small Magellanic Cloud dust.l Drain3l2003 ). but nevertheless it 
is clear that the scattered light observations require much lower 
densities than those implied by [S II] ratios. The two measure- 
ments can be reconciled if the gas is highly clumped, so that 
most of the luminosity is coming from high-density clumps, 
whereas the mass and the scattering cross-section are domi- 
nated by low-density gas. 

While a detailed modeling of all observables is beyond the 
scope of this paper, we use a toy model in which the mass of 
the emitting gas at each density is a power-law function of the 
density, with power-law index -a between Wmin and n,nax, to 
estimate the clumping factor Since the [S II] line ratios are 
usually observed to be between the low-density and the high- 
density asymptotes, values of Wmax = a few times the critical 
density are required; we use n^ax = 10"' cm"-'. At the same 



time, for 1 < a < 2 the minimal density is constrained to be 
Umin = (1 - ct){ne) /o; by the scattering observations. With these 
constraints, the clumping factor is 

{\-af ( nmaxV'" 

K = : ■ (4) 

a{2-a) J 

For example, for a = 1.5 and («(.) = 0.016 cm"^, for each 
lO**' erg s"' of H/3 emission, the mass of the emitting gas is 
Mg = 7 X IQ^Mq. This estimate can only be considered very 
approximate, since the derived mass is quite sensitive to the 
specific assumptions about clumping (for example, it varies by 
2 dex as a varies between 1 and 2). Nevertheless, we point out 
that the standard method of mass determination likely produces 
an underestimate of the true mass and that the scattering obser- 
vations provide a valuable constraint on the physical conditions 
in the NLR. 

In short, we see compelling evidence that the NLR is clumpy. 
As a result, it is difficult to estimate robust gas masses, and thus 
difficult to determine what fraction of the gas may be expelled 
by potential AGN outflows. 



5. TWO CANDIDATE DUAL OBSCURED QUASARS 

Many recent surveys have identified potential dual active 
galaxies (i.e., two active galaxies with ^ kpc separations) as 
narrow-line objects with m ultiple v elocity peaks in the [O III] 
line in SDS S spectra (e.g.. lWang et al. 2009; Li u et al.ll20T0b( 
Smith et al. 2010), as w ell as from the DEEP2 redshift survey 
dGerke et al.. 2007; Comerford et al...2009b.) . Other candidates 
have been identified based on spatially offset nuclei dBarth et alj 
l2008l;IComerford et al.l200"9al) . There are two intriguing objects 
in this sample that may contain dual AGNs. 
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Fig. 9.— r-band image of the dual AGN SDSS J0841+0101. The primary 
object (A) is indicated with the black an'ow, red aiTow points North, and the 
scale bar is 10" long. Separation between the two AGNs is 3.8 " (7.6 kpc). 

The first is SDSS J1356+1026 (Fig. |3]l, which has two clear 
continuum sources, each with associated high-ionization [O III] 
emission. Their separation is ^ 2.5 kpc (l."lV_This object was 
highhghted as a potential dual AGN by both lLiu e t al. (2010b) 
ba sed on multiple velocity peaks in the SDSS spectrum and 
by lFuet al.l(l2010h from Keck AO imaging. We have recently 
shown that ~ 10% of the double-peaked narrow-line c andidates 
also ha ve spatially resolved dual continuum sources ('L iu et all 
I2010ah . It seems natural that two galaxies would contain two 
BHs. On the other hand, there well may be a single radiat- 
ing BH that is illuminating all of the surrounding gas. Unfor- 
tunately, our long-slit spectra do not include [O II] or [S II], 
which would give us a handle on the electron densities, and 
thereby whether a single ionizing source is plausible. Given the 
projected separation of 2.5 kpc, if we assume that there is a sin- 
gle ionizing source associated with one of the two continua, we 
would expect to see the ionization parameter decrease by a fac- 
tor of 6 between the two targets. In fact, the [O III]/H/3 ratios 
are within 10% of each other, as are the [O III] fluxes. On the 
other hand, the very high ionization parameter seems to extend 
over the entire nebulosity 10 kpc). Of course, the accret- 
ing BH may sit between the two continuum sources. Definitive 
proof requires the detection of X-ray or radio cores associated 
with each continuum source. 

SDSS J0841-I-0101 shows much less ambiguous evidence for 
a pair of accreting BHs, with a projected separation of 3 "8 (7.6 
kpc; Fig. |9j. It would not be included in double-peaked sam- 
ples assembled from the SDSS because the separation on the 
sky between the two components is larger than the SDSS 3" 
fibers. Nevertheless, the component separations are com para- 
ble to those in the Liu et al. sample. Liu et alj ( l2010a) show 
that the double-peaked samples are probably dominated by sin- 
gle AGNs. These observations highlight that we are likewise 
missing dual AGNs with slightly larger separations. 

As is apparent from Figure [TO] the two AGNs are strik- 
ingly similar in spectroscopic properties. The [O III] lumi- 
nosities (« 10**^ erg s"') agree within < 0.1 dex, and the 
[O m]/H/3 ratios (-- 10) agree within 5%. The only clear dif- 
ference is in the linewidths. The primary galaxy (A) has a 



FWHM[oin]= 430 km s"', while the companion AGN (B) is 
narrower, with FWHM[oin]= 330 km s"'. This difference most 
likely reflects the fact that A, with a stellar velocity disper- 
sion of (T* = 214 ±29 km s"' is more massive than B, with 
fj* = 150 ±40 km s"'. Taken at face value, this difference in dis- 
persions corresponds to a difference of a factor of nearly 10 in 
BH mass between the two galaxies. Accordingly, if the [O III] 
luminosity tracks the bolometric luminosity, then apparently B 
is accreting 10 times closer to its Eddington limit than A. 

Alternatively, there may be only a single radiating black hole. 
If there is only one quasar - in galaxy A - then we consider two 
scenarios. The first is that the quasar in A is unobscured as seen 
from B, so that the galaxy B is photoionized by the central en- 
gine in A. If we assume that most of the NLR emission in A is 
produced at a distance < 1 kpc from the nucleus, then in order 
to preserve the ionization parameter (as evidenced by the sim- 
ilar spectra of A and B), the difference in electron density be- 
tween the two galaxies would have to be a factor of > 60. While 
the dust particles in galaxy B may scatter quasar spectrum, this 
emission cannot dominate the observed spectrum (otherwise we 
would see a broad-line AGN in source B). The resulting esti- 
mates of the emerging equivalent width of the emission lines 
suggest that this scenario is possible, but has to be quite tuned 
in order to fit observations. The second scenario is that galaxy 
B is located along the obscured direction, just like the observer, 
but scatters some of the A's [O III] emission. However, in this 
scenario the ratio of [Olll] fluxes of B and A corresponds to the 
fraction of photons that B intercepts, ^ (2/7.6)2/47r = 0.006, 
contrary to the observed similarity of fluxes. In conclusion, the 
picture of a single active black hole producing two objects with 
similar fluxes and ionization parameters appears unlikely. 

6. DISCUSSION & SUMMARY 

We are looking for direct signs of feedback in the two- 
dimensional spatial extents and kinematics of the NLRs of a 
sample of luminous obscured active galaxies. Our conclusions 
are mixed. On the one hand, we see clear evidence that the 
AGN is stirring up the galaxy ISM. On the other hand, we do 
not see signs of galaxy-scale winds at high velocities. However, 
as we argue below, perhaps this is unsurprising. 

We see two distinct signatures of a luminous accreting BH 
on the ionized gas in these galaxies. The NLRs are much more 
extended at these high luminosities than in lower-luminosity 
Seyfert galaxies. In fact, the AGNs are effectively photoioniz- 
ing gas throughout the entire galaxy. This alone means that the 
AGN is heating the ISM on galaxy-wide scales. The impact of 
the AGN is more directly seen in the kinematics. We see very 
few ordered radial velocity curves; instead the velocity distribu- 
tions are typically quite flat even at large radius. Perhaps even 
more striking is that the gas velocity dispersions are high out to 
large radius. As we have argued, not only do inactive galaxies 
uniformly show a drop in gas (and stellar) velocity dispersion 
at large radius, but even in ultra-luminous infrared galaxies the 
gas velocity dispersions are observed to drop at large radius. 
We cannot therefore attribute the gas stirring to gravitational 
effects such as mergers. It is most natural to implicate the ac- 
creting BH. 

On the other hand, overall the velocities we observe in the 
NLR gas are not very high (a few hundred km s"'). Taken at 
face value, our crude estimates suggest that very little of the 
ISM is moving fast enough to escape the galaxy, although a 
clumpy NLR complicates our ability to estimate this fraction 
robustly. In only one case do we see the spectacular outflow- 
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Fig. 10. — The binary narrow-line quasar serendipitously discovered in our spectroscopic campaign. Component A was observed by the SDSS and is the more 
massive of the two systems. Left: Stellar velocity dispersion fits using the best weighted fit of an F2, Kl, and K4 star. These fits provide both a measurements of cr* 
and also a continuum subtraction. Flux densities are normalized to the continuum. Right: Fits to the H/3 and [O iii] emission lines using multi-component Gaussians; 
see text for details. 
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ing nebulosity one might imagine in thinking of AGN feedback 
(SDSS J1356+1026). Before we can rule out that any gas is 
unbound from these galaxies, however, we should consider the 
impact of projection effects, potential observational biases, and 
some theoretical expectations. 

Our observations suggest that ionized gas is ubiquitous 
within the galaxy but rare at larger (e.g., 10 kpc) scales. As 
explained above, the observed ratio of obscured to unobscured 
objects leads us to assume an ionization cone opening angle of 
^ 120°. With such a large opening angle, we would expect our 
slit to intercept the NLR nearly all the time, as we observe. On 
the other hand, we see extended gas on 10 kpc scales in only 
one case. Furthermore, the HST continuum images show ex- 
tended emission on these large scales, but with a much smaller 
opening angle of 20-60°. Similarly, we have visually inspected 
the most luminous obscured AGNs from the Reyes et al. sam- 
ple with 0.16 < z < 0.3 and found evidence for small open- 
ing angles from the broad-band images (which have significant 
[O III] light in the r band). Probably we are seeing the effects 
of surface brightness dimming at the outer reaches of the bi- 
cone. Although the true opening angle is large (120°), only 
a much narrower inner cone can be observed at 10 kpc. Tak- 
ing the smaller opening angles, we expect to see extra-galactic 
extended gas only 20-40% of the time. That fraction is not 
inconsistent with the number of objects that we observe with 
emission line regions extending beyond their host galaxies. 

Projection effects also preferentially bias us against detect- 
ing the true outflowing velocities. These are obscured objects, 
and on large scales we see evidence for ionization cones in the 
HST continuum imaging. We thus expect the largest acceler- 
ations to occur in the plane of the sky. We perform a Monte 
Carlo simulation in which the NLR is modeled as a biconical 
outflow with constant velocity as a function of radius, assuming 
different opening angles for the bicone (Fig. fTTT i. We sample 
random lines of sight outside of the bicone, and find that while 
the intrinsic velocity is uniformly high, we only expect to ob- 
serve high (e.g., approaching escape) velocities a small fraction 
of the time. 

These simulations take into account only the bias introduced 
by projection effects and assume constant velocity and uniform 
emissivity within the bi-cone. We also considered more real- 
istic models, in which velocity varies as a function of distance 
(-2 < dlogv/dlogr < 1) from the center, mass conservation is 
satisfied and the emissivity correspondingly declines as «^(r). 
Due to the decline of emissivity in these models, at a projected 
distance d from the center the observed brightness is dominated 
by the location physically closest to the center (that is, r = d), 
and this gas is moving exactly in the plane of the sky exacer- 
bating the projection bias. In the case v(r) cx the emissivity 
may be uniform, but the integral along the line of sight is dom- 
inated by gas that moves slower than the gas at d because of 
the declining velocity profile. Therefore, in these more realis- 
tic situations we find a radial velocity distribution that is more 
peaked at zero than shown in Fig. [TT] Thus, while we observe 
small escaping fractions, once projection effects are accounted 
for, the observations may be consistent with high velocities in a 
large fraction of the gas. 

Finally, there is the possibility that the outflows operate pre- 
dominantly on small scales. In local low-luminosity sources 
outflows are observe d only with in the inner hundreds of pc 
(e.g.,|Crenshaw & Kraemei 120001) . In addition, recent simula- 
tions bv iPebuhr et al.l (12010 ') suggest that BHs do self-regulate 
their own growth but do not generate galaxy-wide outflows. 



Of course, other factors may be at play as well. There is 
the possibility that some fraction of the ioni zing photons have 
escaped the galaxy (e.g.. iNetzer et al.l 12004). or even that we 
are seeing galaxies in some pre-outburst phase, as may be ex- 
pected if obscured accretion tends to accom pany the late stages 
of merging and star f ormation activ itv (e.g..iMihos & HernquistI 
1994; San ders & Mi rabel 1996; Ca naUzo & Stocktonll200lT 

It is interesting to compare with simulations of galaxy-scale 
out flows. We start with the work of P r ogaetal. (2008) (see 
also lProgall2007l; iKurosawa & Progall2009l) . These simulations 
focus on smaller scales than those we probe, extending no fur- 
ther than 10 pc. However, it is at least a starting point for com- 
parison. The simulations include radiative heating by both an 
accretion disk and an X-ray corona, and look at the impact of 
varying the density and temperature structure, as well as rota- 
tion, of the gas. We highlight a few generic conclusions from 
their studies that are very relevant to our work. First of all, 
the final flow includes both an equatorial inflow and a bipolar 
outflow. Consistent with our work, the opening angle of the 
outflowing cone can be quite wide (up to 160°). Also inter- 
esting to note is that the outflows can be dynamic, clumpy as 
we observe, and with multiple temperatures (ranging from the 
10"* K gas observed here all the way to X-ray emitting temper- 
atures). It remains to be seen whether the outflows on pc scales 
will propogate to larger (galaxy-wide scales). 

A recent study by Hopkins et al. (in preparation) of outflows 
driven by AGNs in numerical simulations demonstrates several 
surprising similarities to the kinematics of the ionized gas we 
see in our study. The observations suggest that outflows are 
clumpy because the measurements of rms density and the mean 
density are highly discrepant. The simulations suggest that out- 
flows are clumpy because they are subject to Rayleigh-Taylor 
instabilities. Furthermore, the rate of the decline of mean den- 
sity with distance from the center seen in scattering observa- 
tions is similar to that seen in numerical simulations where the 
motion of the gas becomes ballistic at large distances. The 
masses and the velocities of the outflows that we find are quite 
similar to those seen in numerical simulations, and although the 
kinetic energies of the outflows (^ 10"*^ erg s"') are just a small 
fraction of the total energy output of the AGN, the simulations 
suggest that the wind is in fact driven by a much stronger cou- 
pling of the AGN output to the gas. The small kinetic energies 
that we see at this late (^ 10^ years) stage are simply left-overs 
after much of the energy was efficiently radiated by the outflow. 
While these qualitative similarities are very encouraging, the 
specific mechanism responsible for coupling of the black hole 
output to the gas on much smaller spatial scales (which then de- 
velops into the relic outflow we see now) remains unidentified. 

In short, it is clear that the presence of the AGN at the galaxy 
center impacts the entire galaxy. Whether significant mass out- 
flows are driven, particularly in the radio-quiet regime consid- 
ered here, remains an open question. The next step for this type 
of an alysis is akeady und erway. Integral-field unit observations 
(e.g.. lVillar-Martfn et al fcOlOi) . particularly with a wider wave- 
length coverage, will remove some of the ambiguities we strug- 
gle with. 

We thank G. Novak for numerous interesting discussions, 
and P. Hopkins for sending us a manuscript in advance of pub- 
lication. We thank the referee, Sylvain Veilleux, for a very 
prompt, careful and helpful report that significantly improved 
this manuscript. Research by A.J.B. is supported by NSF grant 
AST-0548198. 
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Fig. 11. — To illustrate the magnitude of the projection effects, we present the results of two Monte Carlo simulations of observed velocity distributions for 
a conical outflow with an opening angle of 60° {solid histogram) and 120° {dotted histogram). The line of sight is constrained to fall outside of the cone but is 
otherwise drawn from the appropriate random probability function. The outflow is assumed to emit uniformly and to have a constant outflow velocity Vacuj, but we 
measure the (smaller) radial velocity Vobs- On the left, we show the distribution of velocities from both bicones, while on the right we show only the approaching 
bicone. Because many of the streamlines of the gas lie close to the plane of the sky, the observed velocities are biased to be significantly smaller than the real ones. 
This figure demonstrates that even if we detect only a small fraction of the gas at high velocities approaching the escape speed, a large fraction of the gas may 
actually be escaping. 
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APPENDIX 

This appendix includes all of the two-dimensional infomation for all galaxies that are spatially resolved in our observations (Figs. 
12-20). Note in particular the high velocities and dispersions at large radius. 



16 







' ' ' I ' ' ' I ' ' ' I 

J0B41+0101 PA=44 



_L 



_L 



_L 





200 



E 



-200 
750 



J 500 



250 




CQ. 10 



£1 C) 



(I U 



H \ 1 1 \ h 



o 



o 



4200 4400 4600 4800 5000 
Rest Wavelength (A) 

10 



o o 

H \ 1 \ 1 \ 1 h 



GREENE, ET AL. 



0.5 







T 



3 
- 



500 



250 



10 



_L 



_L 



_L 



200 — I — I — ^ 



100 



-100 



-200 1 1 h 



o o o 



-\ — I — I — h 



_L 



4200 4400 4600 4800 5000 
Rest Wavelength (A) 




T <j' 4 
° t " T 

o o o ° ° 

H — I — \ — I — I — I — I — I — I — ^ 



Radius (arcsec) 



Radius (arcsec) 



Fig. 12. — Resolved spectroscopic measurements for SDSS J0841+0101. In the top panel we show the central spectrum, extracted within the inner ~ 1" arbitrarily 
normalized. The next panel shows the flux in the [O iii] line (filled black squares) and in the continuum (open stars) as a function of radius (10"'^ erg cm"- s"'). The 
velocity panel shows the velocity offset between [O iii] and the "systemic" velocity as measured from H/3 in the nuclear spectrum (filled black circles), the velocity 
dispersion in the gas (open black circles), and the maximum velocity (measured at 20% of the line maximum) to the red (red open circle with cross) and to the blue 
(blue open circle with dot). Finally, the ratio of [O !n]/H/3 is plotted (black crosses). A ratio of [0 iii]/H/3 of three is noted with the dashed line. 
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Fig. 14.— Resolved spectroscopic measurements for SDSS Jl 124+0456 and SDSS Jl 142+1027. Symbols as above. 
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Fig. 15. — Resolved spectroscopic measurements for SDSS J1253-()341 and SDSS J1413-0142. Symbols as above. 
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Fig. 17. — Resolved spectroscopic measurements for SDSS J2056+0057. Symbols as above, shown for each slit position. 
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Fig. 18. — Resolved spectroscopic measurements for SDSS J2126+0035. Symbols as above, shown for each slit position. 
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Fig. 19.— Resolved spectroscopic measurements for SDSS J1507+0029 and SDSS J2212-0944. Symbols as above. 
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